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ABSTRACT 

Alfven turbulence caused by statistically isotropic and homogeneous primordial magnetic field induces 
correlations in the cosmic microwave background anisotropies. The correlations are specifically between 
spherical harmonic modes ai-\^ m and a^+i jm . In this paper we approach this issue from phase analysis of 
the CMB maps derived from the WMAP data sets. Using circular statistics and return phase mapping we 
examine phase correlation of A£ = 2 for the primordial non-Gaussianity caused by the Alfven turbulence 
at the epoch of recombination. Our analyses show that such specific features from the power-law Alfven 
turbulence do not contribute significantly in the phases of the maps and could not be a source of 
primordial non-Gaussianity of the CMB. 



1. INTRODUCTION 

After the release of the 1-year WMAP data for anal- 
ysis (Bennett et al. 2003a, b,c; Hinshaw et al. 2003a, b), 
one of the most exciting area of investigation is in the 
statistical characterization of the CMB signal. After the 
WMAP science team's report on the observational con- 
straint on the quadratic non-Gaussianity of the anisotropy 
of the CMB (Komatsu et al. 2003), Chiang, Nasel- 
sky, Verkhodanov & Way (2003); Coles, Dineen, Earl & 
Wright (2004); Naselsky, Doroshkevich & Verkhodanov 
(2003, 2004); Vielva et al. (2004); Park (2004); Erik- 
sen et al. (2004); Hansen et al. (2004); Larson & Wan- 
delt (2004) have reported detections of non-Gaussianity 
by various methods for the CMB maps derived from the 
WMAP data. Unfortunately, the origin of such non- 
Gaussian features, detected by different methods, is still 
unknown. In Chiang, Naselsky, Verkhodanov & Way 
(2003); Naselsky, Doroshkevich & Verkhodanov (2003, 
2004); Eriksen et al. (2004), they point out that the 
non-Gaussian features might come from the foregrounds, 
whereas Hansen et al. (2004) argue in favor of systematic 
effects. Moreover, different methods detect various proper- 
ties of non-Gaussian features from localized peculiar spots 
(Vielva et al. 2004) to global north-south asymmetry in 
the WMAP signal and COBE map as well (Eriksen et al. 
2004). 

There is another point of view on the issue of the 
WMAP non-Gaussianity. It comes from the theory of tur- 
bulence in magnetohydrodynamics (MHD). Using the so- 
called extended self-similarity (ESS) method, Bershadskii 
& Sreenivasan (2003, 2004) point out that the statisti- 
cal properties of angular increments ST r = T(R + r) — 
T(R), where r is the vector connecting two pixels R + r 
and R of the map, have the following relation in MHD 
(\5Tr \ p ) ~ r^p for high order moments p = 4, 6 ... 12. For 
those moments, the ESS is clearly detected both in the 
COBE and the WMAP data : (\8T r \ p ) ~ (<5T r 4 )«» where 



Cp — § + 1 — (f) 3 ■ Moreover, in Branderburg, Enqvist & 
Olesen (1999); Durrer, Kahniashvili & Yates (1998); Mack, 
Kahniashvili & Kosowsky (2002); Subramanian, Seshadri 
& Barrow (2002); Chen et al. (2004) it was shown that pri- 
mordial magnetic field, one of the most significant relics 
from inflation, can generate vorticity of Alfven waves be- 
fore and during the epoch of hydrogen recombination. It 
interacts with the CMB photons, producing non-adiabatic 
tail of the CMB anisotropy and polarization. Assuming 
statistical homogeneity and isotropy on primordial mag- 
netic field, it is shown that such Alfven turbulence shall 
have non-Gaussian properties, because of quadratic depen- 
dence of the vorticity amplitude on the magnetic strength 
B. Namely, such a magnetic field induces correlation be- 
tween the ae,-x, m an d a>l+i,m multipole coefficients of the 
CMB temperature anisotropy expansion by the spherical 
harmonics. 

An intriguing issue then follows: has WMAP observed 
the cosmological Alfven turbulence? Recently this issue is 
discussed in Chen et al. (2004), which exploits the corre- 
lations between a^_i jm and a^_|_i, m of the CMB caused by 
the vorticity of Alfven waves in order to put constraints 
on the strength of the magnetic field: |B| < 15 nG for the 
spectral power index n = — 5, and |B| < 1.7 nG for the 
spectral power index n = — 7. 

In this paper we discuss this issue from another aspect: 
could the detected non-Gaussianity from the WMAP data 
be related to the Alfven turbulence? The basic idea is that 
even if the Alfven turbulence is small in amplitude in the 
CMB signal, it shall manifest itself in the WMAP maps 
between the phases 4>£-i jm and <^+i, m - Phase correlation 
between Fourier modes has been investigated in relation 
to large-scale structure formation of the universe (Scher- 
rer, Melott & Shandarin 1991; Chiang & Coles 2000; Chi- 
ang, Coles & Naselsky 2002) and also applied as a test on 
non-Gaussianity (Chiang, Naselsky & Coles 2002; Chiang, 
Naselsky, Verkhodanov & Way 2003; Coles, Dineen, Earl 
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& Wright 2004) based on the random phase hypothesis as 
a practical definition of Gaussian random fields (Bardeen 
et al. 1986; Bond & Efstathiou 1987). Naselsky, Novikov 
& Silk (2002) use neighboring phase correlations to extract 
extragalactic point sources. 

Assuming that the CMB signal is composed of pure 
Gaussian signal and sub-dominant vorticity tail, we use 
circular statistics on phases for such a signal in order to 
place constraint on the power spectrum of Alfven turbu- 
lence at each multipole number t. We will show that cir- 
cular statistics allow us to decrease contamination of the 
dominant (Gaussian) part of the signal in order to inves- 
tigate possible contamination from the non-Gaussian sub- 
dominant part. 

2. CIRCULAR STATISTICS OF THE ALFVEN TURBULENCE 

For statistical characterization of temperature fluctua- 
tions on a sphere we express each signal (either CMB or 
foreground components) as a sum over spherical harmon- 
ics: 

oo I 

AT(e, v ) = J2 E WmW^v^e^), (i) 

£=0 m=-t 

where \ae m \ and (j>e m are the moduli and phases of the 
coefficients of the expansion. 

Homogeneous and isotropic CMB Gaussian random 
fields (GRFs), as a result of the simplest inflation 
paradigm, possess Fourier modes whose real and imagi- 
nary parts are Gaussian and mutually independent. The 
statistical properties are then completely specified by its 
angular power spectrum C| mfc , 

(oMovcr 4 ^^'- (2) 

In other words, from the Central Limit Theorem their 
phases 



.j.cmb 
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are randomly and uniformly distributed at the range 
[0, 2n]. We hereafter denote with G as the pure Gaussian 
tail of the CMB signal. For the combined signal 

a-em — Gem + Vim, (4) 
where Gi m is the Gaussian tail and Ve m is the vortex tail. 
We can write down 

ai m = \at m \ exp(i*£ m ), 

where 



(3) 



\ai n 



tan $>e m = 



\Glm\ Z + \Vlm\ Z + 2\G tm \\Vt m \cOB(*f m - * , nl 

| Gj m | sin ^f m + | V tm | sin 
\G lm \ cos *f m + \Vem\ cos *Y m 



(5) 

v ); 

(6) 



and &f m , ^Jrn are the phases of Gaussian and vortex at 
the (£, to) harmonics, respectively. 

In Durrer, Kahniashvili & Yates (1998) (hereafter DKY) 
it is shown that for the CMB anisotropy produced by the 
Alfven turbulence, the properties of the power spectrum 
and phases are different from the adiabatic modes, which 
corresponds to the Gaussian tail of the signal. Namely, 
for the power spectrum CV(to) = (o^ m a| m ) and the auto- 



correlator Dt{m) = (a 
C/(m) = A 



l,m a e+l,mi 
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from DKY we get 



D t {m) 



2£ 4 + 4l 3 - l 2 - 31 + (6 - 2£ - 2£ 2 )m 2 
{2t-\){2t + Z) 
2A {T{-n- 1)T(£ + 5 + 3 ^ 



| n + i| r(-f) 2 r 



2' V 2*2' 

{£ + m + l){£ - m + l){£ + m){£ - m) 



2 ' {£-!)(£ + 2) x 
(7) 



{21- l)(2£+l) 2 (2£ + 3) 

where A is the normalization constant and —7<n< —1 
is the power spectral index of the Alfven turbulence. Be- 
low we discuss the models with n = —3, —5 and —7, for 
which the mean power spectrum over m direction in Eq. (7) 
are C(, De oc £ n+3 . Note that the model of n = —5 
corresponds to the power spectrum Ce oc £~ 2 , which at 
the multipole range of 5 < £ < 20 has the same feature 
as that of the CMB for adiabatic perturbation. For the 
model of n = — 7, on the other hand, the power of ver- 
tex tail is Ce oc and increase rapidly if I — > 0, which 
is typical for the diffuse extragalactic foregrounds. For 
n = — 3 the vertex power Cg <~ const mimic the power of 
extragalactic point source component. Therefore, simple 
statistics based on the estimation of the D^-th moment 
from the WMAP data is potentially misleading, because 
of possible contribution from such kind of foregrounds as 
sub-dominant components to the signal. 

The method we propose in this paper is based on the 
statistical characterization: non-Gaussianity of the vortex 
part of the CMB signal and related with correlation of the 
phases of (I — 1, to), (I + 1, to) harmonics. The correlation 
vanishes for pure Gaussian CMB signal, whereas it is sig- 
nificant for the vortex component. We therefore, introduce 
some specific functions of phases which minimize the con- 
tribution from non-correlated Gaussian signal and maxi- 
mize the non-Gaussian tail of the phases. For this purpose 
we use trigonometric moment statistics to counter the cir- 
cular nature of phases (Fisher (1993), see also Naselsky, 
Doroshkevich & Verkhodanov (2003)). Let us define the 
following trigonometric moments: 



1 



e-i 



cos(^_i, m - Ve+i 



S(£) = J^l J2 sin(*«-i,n» - */+i, m ); 

m— 1 



r 2 (£) = C 2 (£) + S 2 {£); 



R = 



'-max *-min 



+ 



(8) 



The reason for such statistics is clear. Simple algebra leads 
to the following properties of the phases: 

COs(*£_i, m - *M-l, m ) = 

\gi-i, m \\ge+i,m\ cos(*f_ l TO - *f +1>m ) 

+ h-l,m|h+l,m| COs(*f_ liTO - ^Y+l,m) 

+\9e-i, m \\ve+i, m \ cos(*f_ 1)m - ^Y+i.m) 



+ |5m,m||^-l,m| COs(*f +1 ^ 



t-\,m) 



(9) 



where \g tm \ = \Ge m \/\a im \,\v em \ = \V em \/\a tm \. As one 
can see from Eq.(9) the first term is proportional to 
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to) = (a l _ 1 m a} +1 m ) moment for the pure Gaussian 
tail of the signal, which should be vanish, the second term 
is D^{m) for the vortex, and last two terms correspond 
to correlations between G and V, which should be statis- 
tically negligible after summation over to. If the vortex 
component is sub-dominant, then in Eq.(6) and Eq.(9) 
\Ve m \ < \Gt m \ for all £,m and \a im \ is given by Eq.(6). 
This means that \gt m \ — 1 and \vi m \ <C 1. However, 
because of finite number of to modes and especially for 
low £ range, correlations can manifest themselves between 
(£ — l,m), (I + l,m) harmonics. In order of magnitude 
this effect can be estimated by the following way. For 
pure Gaussian signal with non-correlated phases 4>t m one 
obtains 

i i 
1 x— ^ 1 1 . , 1 

-^cos<^ ro ~-=; -^sin^ m ~-=. (10) 

m— 1 v m— 1 v 

Then, for the moduli \ag m \ from Eq.(6) one gets 

Km| 2 — |G£ m | 2 + \Vi m \ 2 

if \\V im \^\G im \/VI; 
\a £m \ 2 ~ \G lm \ 2 + 2\G tm \\V lm \ cos(*f m - *L), 

if « |G /m |/V£ (11) 

For the asymptotic \Vt m \ » \Gi m \/\/t, 

1 1 1^-1, m||^+l,m| 

Vf £i ^ |G^_i, m ||G^+i, m | 



C(£) 



+ 



^ |G*-i, m ||G /+1 , 



■COS^-Lm "*/+!, J- (12) 



For the correlated phases of the vortex component, the 
last term in Eq.(12) is larger than the second and, if 



*l-l,m ~ then 



(13) 



where is the power spectrum of the Gaussian tail. 
For highly correlated vortex perturbations, therefore, 

the contribution to the C function is significant if ^fv^ > 

G p 



1. 



At £ = 100, Mi > 0.1, which is exactly the DKY 

G, 



criterion for estimation of the magnetic field amplitude 
for different values of the spectral index n. However, if 
the observable data sets covering the range of multipoles 
up to the Silk damping scale £d for the vortex perturba- 
tions, £d ~ 500 (Durrer, Kahniashvili & Yates 1998), then 

non-Gaussian features could be detectable for =k > 0.045. 

G\ ~ 

Moreover, one interesting feature comes from power spec- 
tral index dependence of the vortex perturbations men- 

tioned above. For n = —7, Dg cx £~ 4 , while G e has the 
same asymptotic £~ 2 as the ACDM WMAP best fit model 
at 5 < £ < 30. For that model of the power index, the 
most stringent constraint should come from statistics of 
the low multipoles from the WMAP data, then from high 
multipole range £ ~ £ d and simple estimator could be 6 -f- 7 
times bigger than in DKY: B% ~ 0.6 + 0.7. 

G e 

We would like to point out that the above-mentioned 
properties of the C statistics provide natural explanations 



in terms of the cosmic variance limit of error bars for any 
CMB experiments. Assuming no systematic effect present 
in the data, for low multipoles £ < 100 the cosmic vari- 
ance limit corresponds to 5Ci/Ct ~ (/ s ky^)~ 5 , where / s k y 
is the sky coverage of the observation. The whole sky 
WMAP coverage gives / s k y = 1 . If some part of the CMB 
power spectrum is related to vortex perturbations, then 
(Chen et al. 2004) 

2 



r(-a±i) 



D t . (14) 



Ce = G e +C e ~G e + \n + l\ r( _„ } 

2 2 

Using the definition 8Ce/Cg — (Ci — G e )/G e , which cor- 
responds to the contribution of the vortex perturbations 
to the power spectrum, we have the following threshold of 
detectability 



II 



r(-^±i 



D t /G t > -= 



1 



(15) 



r(-f) 

which is in perfect agreement with our estimation Eq.(13). 

Let us briefly discuss the properties of the S statistics. 
Similar to Eq.(9) one can prove that sin(\I/£_i jTO — \I^+i )Tn ) 
is given by Eq.(9) with transition from cosine to sine func- 
tion. However, because of sin dependency, if phases of the 
non-Gaussian tail are highly correlated, we have degra- 
dation of the \ve-i, m \\ve+i, m \sm(®Y_ ltm - $Y+i, m ) term, 
even if \Ve m \ » \Ge m \/y/£- So, if the phases of the vor- 
tex perturbations ^Y-i m ~ ^Y+i m) are highly correlated, 

for the S statistics we get S <~ I /VI. Thus, the pres- 
ence of vertex perturbation (as all non-Gaussian signals) 
in the CMB data manifest itself as asymmetry of the C-S 
statistics, while for the pure Gaussian signals it should be 
symmetrical one. In order of magnitude this asymmetry 
for each mode £ is 

S(£) 1 



A{£) = 



C(£) 



G, 



< 0.5, 



(16) 



for highly correlated phases of vortex perturbations. Ad- 
dition to Eq.(16) we can define global asymmetry of the 
signal 



A 



(17) 



£/C(*) ' 

In terms of circular statistical variables (Fisher 1993), 
this global asymmetry corresponds to the mean angle 
6 = tan -1 A g for orientation of all harmonics (£, m) in 
the CMB map. Needles to say for vertex perturbations 
with n — — 5 global asymmetry is expected to be in order 
of magnitude 

and can be detectable, if =#4nax 3> 1- At the end of 

G t 

this section we would like to mention one peculiar asymp- 
totic of the A{£) and A g parameters related with correla- 
tions of phases. From Eq.(12) one can find that correla- 
tion &Y-i,m - *£+i.m = 7r/2,37r/2 leads to A(£),A g > 1. 
So any peaks of the function A(£),A g » 1 are the mark 
of 7r/2,37r/2 correlations for the non-Gaussian tail of the 
CMB signal. Below in addition to definition Eq.(8) we 
introduce the phase difference 



D em (A£) = tt £ 



(19) 
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for even At and 

D im (At) = y e +Ae,m - ^e-At-i, m (20) 
for odd A£ and corresponding circular moments Sa^(^) 
and Ca^(^) similar to Eq.(8). These variables allow us to 
describe properties of the CMB phases for any A£ separa- 
tions of the modes. 

3. CIRCULAR STATISTICS OF THE ILC, FCM AND WFCM 
MAPS DERIVED FROM THE WMAP DATA 

In this section we apply the S, C statistics for the three 
high resolution whole-sky maps: the Internal Linear Com- 
bination Map (ILC) by the WMAP science team, the TOH 
foreground-cleaned map (hereafter FCM) and the TOH 
Wiener-filtered map (hereafter WFCM). The ILC and the 
FCM are obtained from a weighted combination of the 
WMAP 5 maps at frequency 23, 33, 41, 61 and 94 GHz in 
order to separate the microwave foreground. The WFCM 
is the map after foreground residue cleaning by Wiener 
filtering. All the maps formally claim resolution up to 
£ = 1024, but none of them is correct for investigation of 
the CMB properties. The ILC map, as is mentioned at the 
WMAP website, is applicable for the foregrounds property 
investigation at the range of multipoles £ < 300. The FCM 
and WFCM are smoothed for the multipole range I > 300 
by a Gaussian window function and the CMB signal is 
completely erased by the smoothing. However, looking at 
the V and Q bands of the WMAP data and combining 
all the foregrounds (synchrotron, free-free, dust emission) 
and the point sources from the WMAP catalog, one can 
find that this map reproduces the CMB signal outside the 
Kp2 Galactic cut mask extremely well at the range of mul- 
tipoles £ < 50. 

The FCM and WFCM have well-defined scale of smooth- 
ing in the image domain. However, smoothing does not 
change at all the phases of the signal even down to the 
limit £ = 1024. So, using FCM and WFCM we can ex- 
amine the phases at the whole range of multipoles. One 
additional reason for using the FCM and WFCM maps 
is that they are tested by ESS method, from which the 
manifestation of the Alfvcn turbulence properties could 
be clearly seen (Bershadskii & Sreenivasan 2003, 2004). 

In Fig.l we plot for the three maps (ILC, FCM and 
WFCM) the S-C phase diagrams for A£ = (£+l)-{t-l) = 
2. For comparison we include the statistics of Gaussian 
random signal, which reproduce our expectation of the 
phase properties. From the second row in Fig.l one can 
see that ILC phases perfectly reproduce the properties of 
Gaussian signal phases in terms of S(^), C(^) and A{£) 
variables. One can see that on the third and the fourth 
row, which corresponds to the FCM and WFCM maps, 
starting from £ <~ 100 the S-C symmetry is broken, and 
while the S follows Gaussian statistics, its C and A{£) 
appear non-Gaussian. 

According to Fig. 1 , the phases of the FCM and WFCM 
are in agreement with the expectation, described in previ- 
ous Section. We have significant increase of the C statis- 
tics, while S statistics are nearly the same as for the Gaus- 
sian signal. However, next three Figures display clearly 
non-Gaussianity. Fig. 2 are the circular statistics S(t), 
C(£) and A(£) for A£ = 1, and Fig.4 and Fig.3 are for 
At = 3 and At = 4 respectively. According to DKY, for 
all At ^ 2 the phases of the signal should not display any 



correlated features. One can also find remarkable symme- 
try in the FCM and WFCM between even and odd At 
statistics (e.g. the C in Fig.l and Fig.4 and the S in Fig. 2 
and Fig.3). 
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Fig. 1. — Circular statistics on phase correlations with A£ = 2. 
The columns (from left to the right) correspond to S(£), C(£) and 
\ogA(£), respectively. The first row is for a Gaussian realization 
(random phases). The second row is the ILC map, the third the 
FCM map and the last the WFCM map. 



signC^i 



-sign C 2 cos( — A£), for even At 



signSAi — sign Si sin(— Ai?), for odd At. (21) 

This symmetry reflects the symmetry of the phases 
DAe(t), namely 



^£+^,m 

for even A£ and 



7T 



At + TT 



^£+A£,m - ^e-A£-l,m = ^"A£ 

for odd At and for all values of m. 



(22) 
(23) 



4. NON-LINEAR STATISTICS 

In this section in order to investigate the properties of 
non-Gaussianity, we apply non- linear statistics r 2 (t) for 
the ILC, FCM and WFCM phases. From Eq.(8) one ob- 
tains 

cos[D /m (A*) - D en (At)] 

— - -+- y ■> 



r\t) 



1 r + EE 



{t-if 



(24) 



That means, that r-statistics are sensitive to m depen- 
dence of the phase difference for given values £ and At. 
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Fig. 2. — The same as in Fig.l, but for At = 1. 



Fig. 4. — The same as Fig.l, but for At = 4. 
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Fig. 3.— The same as in Fig.l, but for At = 3. 



This statistic has wcll-dcfincd properties. If the signal is 
highly correlated and D tm (A£) — > then r 2 (£) — > 1. If the 
phases are non-correlated, however, r 2 {£) — > In Fig. 5 
and 6 we plot r 2 (£) versus I for the same model and the 
same order as in Fig.l. For all even A£ the correspond- 
ing r 2 (£) are similar to Fig. 5 and for all odd A£ the r 2 
statistics follow Fig. 6. 

Like S and C statistics, r 2 (£) statistics reflect non- 
Gaussian features for A£ = 1,2,3..., which allow us to 
conclude that non-Gaussianity of FCM and WFCM does 
not relate with the vortex turbulence at the epoch of re- 
combination. Moreover, in the next Section we will show, 
how the symmetry of the S and C statistics allow us to 
detect the source of non-Gaussianity of FCM and WFCM. 



5. LOCAL DEFECTS OF THE ILC, FCM AND WFCM MAPS 

In order to detect the nature of the non-Gaussianity of 
the maps with different powers of signals, we will use the 
power filter P(£) = l/\at m \ : proposed by Gorski (1997) 
and Novikov et al. (2001). This linear filter transforms 
any maps with ae m coefficients of spherical harmonic ex- 
pansions to maps with the same amplitude (= 1) at each 
mode £, m, but preserving all the phases and their 
correlations. 



N tm = P{£)a im = exp(i* £m ) (25) 
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tipoles up to £ = 1024 in order to show how power filter 
allows us to detect some of the local defects of the maps 
close to the Galactic plane. In particular, one can eas- 
ily see that there are defects in the ILC map which was 
originally obscured by the amplitudes \ai m \. The clear 
cut features is obviously not of astrophysical origin. One 
crude guess for such peculiar shape in the middle of the 
ILC map is that the WMAP science team might have used 
a direct mask at the center of the map where there are 
pronounced contaminations from foregrounds and then 
used Wiener filtering, combining with the best-fit power 
spectrum, to recover the signal in that area. However, the 
morphology is related to phases (Chiang 2001) and such 
mask shall alter the morphology, which manifests itself in 
the phases. On the other hand, one can see foreground 
residues in both the FCM and WFCM maps. Qualitatively 
speaking, these are clear signatures of noii-Caussianity 
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Fig. 5. — r 2 (l) statistics for Al = 2 for (from top to bottom) a 
Gaussian random realization, ILC, FCM and WFCM, respectively. 
Solid line corresponds to r 2 (£) = l~ 1 asymptotic. 





0.0100 

0.0010 
0.0001 



800 1000 



1 .0000 
0.1000 
0.0100 
0.0010 
0.0001 



800 1000 



1 .0000 
0.1000 
0.0100 
0.0010 
0.0001 




Fig. 6. — The same as in Fig. 5, but for Al = 1. 

In Fig. 7 we plot the power filtered ILC, FCM and 
WFCM maps after power filtration for the range of mul- 



Fig. 7.— The ILC (top), FCM (middle) and WFCM (bottom) 
power filtered maps. These whiten images are to display how fore- 
ground cleaning and the residues manifest themselves in phases, 
which are otherwise obscured by the amplitudes a^ m | . One can see 
the clear cut feature in the middle of the ILC map and galactic and 
point source contaminations for FCM and WFCM. Qualitatively, 
these are clear signatures of non-Gaussianity. 



Using Eq.(25) 



one can easily visualize the 



t— ^ .m 



f+ M. }m cros-correlations using simple definition 



G irn {A£) = N e _^ m N; + ^^ m = exp [iD £m (A£)] (26) 
For these maps the common features, like a small in 
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size clusters, seems to be typical (corresponding statistics 
one can find in Naselsky, Doroshkevich & Vcrkhodanov 
(2004)). In order to show how local defects of the maps can 
produce corresponding features, let us introduce the fol- 
lowing model of sub-dominant non-Gaussian signal. Sup- 
pose that we have a set of points (pixels) with coordinates 
9j , 4>j in which the signal looks like a combination of S- 
f unctions 

AT loc (6,<t>)='22A j 5(<l>-(l> j )5(coB0-cas0 j ), (27) 

3 

where Aj are the amplitudes of defects. In order to obtain 
the corresponding a m coefficients, we convolve AT from 
Eq.(27) with conjugated spherical harmonics 

J d(cos6) J d<f>AT loc (6,<j>)Y; m (0,<j>) 
= Y / A i Y; m {0 i ,4u) (28) 

i 

5.1. Peculiarities of Al- even statistics 

To obtain the correlators for S and C statistics, we need 
to know cos[De m (Al)] and sm[De m (A£)], which can be 
found from Eq.(28) as 

C tm {A£) = |c £ _A£ im ||c £+ A£ m |cos [D lm {M)\ 

= (c e _^L m c* e _^ m + c*^^ m c^^ Jj /2 
S lm (Al) = \c t _M im \\c t+ M tm \sm[D em (A£)] 

= (c e _^ mC ;_^ m - c*_^ m c^^ rn ^ jii 

(29) 

Taking into account Eq.(28), from Eq.(29) we have 

to) cos [m(fa - fa)] 

j,k 

x P™^(cos^)P™ ¥ (cos^) 
S em (A£) = ^2,B jk {£,m) sh^rn^ - fa)] 

x P™ ^ (cos 8 3 )P™ ^ (cos 9 k ) 

(30) 

where P™(cos#) are the associated Legendre polynomials 
and 



which 



B jk (£,m) 
'TO 



AjA k 



Al- m+1 ) T (£+M- m+1 ) 



r( 



(31) 



M + to + l)T(£ + M + m + 1) 

Following Naselsky, Doroshkevich & Vcrkhodanov 
(2003), let us describe an ideal situation when 6j = 6 k = 
it/2, but <f>j 7^ fa. For this model in Eq.(30) we get 

l(0) = 



r l 




As one can see from Eq.(32) and Eq.(30) the major part 
of the summation over to is related with to = £ — A! ) for 



PT-^) P T + ^)\m=t-^ = 



o2f-A£_ r (2 
Z 7T 



r(i + Ai) 



cos(—-) (33) 



Thus, the sign of the Ci m (A£) for even Al is determined 
by the cosine terms. Moreover, if Al = 2n, n = 1,2,3..., 
the contribution of the local defects with 9j — it/2 to C 
statistics manifest itself as the changes of the sign, men- 
tioned in Eq.(21). One may ask, why are S and C statistics 
so different, if they depends on 0j in the same way? Taking 
into account the <j>j dependency of Ce m (A£) and Si m {A£) 
one can see that for j = k the cosine terms in Eq.(30) 
(cos [m(fa — fa)]) goes to unity while the sine terms for 
Se, m (A£) is equal to zero. All j ^ k modes in Eq.(30) are 
represented by highly oscillated functions and look like a 
noise in the phase diagrams (see Fig. 1-4). 

5.2. Peculiarities of odd A£ statistics 

Let us discuss the model Al = 2n + 1, n = 1,2,3.... 
There are two possibilities to understand the properties 
of the phase correlations for odd Al. Firstly we discuss 
the Galactic plane sources contamination. In such a case 
similar to Eq.(33) one can find 



i^A/(0)iy(0)L=/- 



A£ 



(ttA£\ 
C0S {—) 



(34) 



and odd Al does not contribute to S and C statistics, if 
0j = it/ 2 for all j. However, if for some of the local defects 
0j ^ it/2, but cos#j <C 1, there should produce significant 
peculiarities in S statistics by the following way. Close to 
the Galactic plane we can expand the Legendre polynomi- 
als using a Taylor series 



Pt-Aii^ej) = pp_ Ae (o) + 



dPV 



-AC 



{cos6 j ) i 



d(cosdj) 
ld 2 Pp_ Ai {cos9 J ) 
2 d(cos0j) 2 



|cos# 3 =0COStfj 



cos<?, =0 COS 



(35) 

Simple algebra allows us to conclude that the dependence 
S^ TO (A^), Ce m (A£) oc cos 2 0j sin(7rA^/2) comes from in- 
terference between the first and the last terms in Eq.(35), 
after their substitution to Eq.(30). Thus, this dependency 
over Al takesplace for some of the multipoles, for which 
I6j < 1, while for 16 j > 1 S tm {A£), C tm (Al) statis- 
tics goes to zero for odd Al. This dependency of the 
S^ TO (A^), Cz m (Al) over 51 can be clearly seen in Fig. 2 and 
Fig. 3, where symmetry S and C statistics are restored for 
high I. However, Se m (Al) as C( m (A£) depends on Al 
in the same way, while in Fig. 2 and Fig. 3 one can see 
significant asymmetry in the S and C statistics. The rea- 
son for that could be localization of the the defects in 
(j) direction close to the Galactic center with the width 
— it/2 <C 5(j)j -C 7r/2 in the Galactic coordinates. 

6. CONCLUSION 

In this paper we examine the primordial non-Gaussian 
feature from Alvan turbulence caused by primordial mag- 
netic field. The distortion of the blackbody CMB spec- 
trum by such homogeneous magnetic field is not discussed 
in this paper. Due to the vector nature of the magnetic 
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field, off-diagonal correlations between spherical harmonic 
modes a^_i iT7J and ae + i, m are induced. To see the At = 2 
correlations, we apply circular statistics of phases to an- 
alyze such non-Gaussian component from the CMB sig- 
nal. We have applied the statistics on the ILC, FCM and 
WFCM maps, using all the avaliable phases. The phase 
information can help us to test some of the properties of 
the signal even when the amplitude of the power spec- 
tra were smoothed by Gaussian window functions starting 
from £ ~ 50 - 200. 

The FCM and WFCM maps have non-Gaussian signa- 
ture registered in the phases, as we have found phase cor- 
relations not only A£ — 2, but also A£ = 1, 3 and 4. Such 
strong correlations are related to the residues from com- 
ponent separation in the Galactic plane and some point 
sources. None of the maps shows specific features from the 
vortex perturbations at the last scattering surface. More- 
over, the detected non-Gaussianity is completely different 
from the vortex perturbations properties. Roughly speak- 
ing, all the modes with At = 1, 2, 3 and 4 show pro- 
nounced correlations, which is not the specific features for 



the vortex. The correlations such as A£ — 4, 8 . . . might 
have been the higher order correlation from A£ = 2, but 
the At = 1 and others indicate that such strong correla- 
tions are related to the residues from foreground compo- 
nent separation. We conclude that the Alfven turbulence 
does not contribute significantly in the phase informa- 
tion. The method we propose for phase analysis is useful 
for the upcoming PLANCK mission, especially for testing 
the properties of foreground cleaning methods. Our final 
remark is that the MHD structure claimed by Bcrshad- 
skii (2003) and Bershadskii & Sreenivasan (2004) in the 
WFCM map is related with Galactic plane residues but 
not with the Alfven turbulence contamination. 



We thank Tegmark et al. for providing their processed 
maps. We acknowledge the use of Healpix 3 package 
(Gorski, Hivon & Wandelt 1999) to produce ai m from the 
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vich et al. 2003) for data analyses and the whole-sky fig- 
ures. 
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